Mode coupling and conversion at anticrossings treated via 
stationary perturbation technique 

Dzmitry M. Shyrokj*] 

Department of Theoretical Physics, Belarusian State University, 
Fr. Skaryna avenue 4, Mensk 220080, Belarus 

Abstract 

Intermodal interactions displayed through the phenomena of mode coupUng and conversion in 
optical systems are treated by means of the Lindstedt-Poincare perturbation method of strained 
parameters more widely known in classical quantum mechanics and quantum chemistry as the sta- 
tionary perturbation technique. The focus here is on the mode conversion at the points of virtual 
phase matching (otherwise called anticrossings or avoided crossings) associated with the maximum 
conversion efficiency. The method is shown to provide a convenient tool to deal with intermodal 
interactions at anticrossings — interactions induced by any kind of perturbation in dielectric index 
profile of the waveguide, embracing optical inhomogeneity, magnetization of arbitrary orientation, 
and nonlinearity. Closed-form analytic expressions are derived for the minimum value of mode 
mismatch and for the length of complete mode conversion (the coupling length, or the beat length) 
in generic waveguiding systems exhibiting anticrossings. Demonstrating the effectiveness of the 
method, these general expressions are further applied to the case of TE„ ^ TM^ mode conversion 
in (i) a multilayer gyrotropic waveguide under piecewise-constant, arbitrarily oriented magneti- 
zation, and (ii) an optically-inhomogeneous planar dielectric waveguide — an example which the 
standard coupled-mode theory fails to describe. 
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I. INTRODUCTION 



Mode coupling and conversion are commonly known and primarily important phenomena 
in fiber and integrated optics, either hindering or fueling the operation of numerous devices 
and elements, e.g., electro-optical and acousto-optical switches and modulators, waveguide 
couples, power splitters, wavelength filters, and others jl| . When a mode of certain configu- 
ration transverses the structure and converts into some other mode, say, due to the externally 
induced perturbation in optical properties or geometrical configuration of the structure or 
due to the imperfectness of the materials used, or when the exchange of energy between the 
guided modes of adjacent waveguides occurs — anyway we face the phenomena of mode 
coupling and conversion. The crucial regime for the mode coupling is in the vicinity of the 
points of virtual phase matching (otherwise referred to as "anticrossings" or "avoided cross- 
ings" — a terminology brought about from quantum mechanics and solid-state physics). 
At anticrossings, the complete mode conversion can be achieved, which makes tailoring this 
regime particularly important for the applications. For a positive example, let me mention a 
promising method for controllable dispersion compensation in photonic bandgap fibers based 
on weak interactions at anticrossings between the core-guided mode and a mode localized 
in an intentionally introduced defect of the crystal 2]. An opposing example of technology 
strains due to the undesirable intermodal interactions is radiation losses caused by the cou- 
pling of the guided mode to radiation modes in a two-dimensional photonic crystal etched 
into a planar waveguide j^l . 

A dominating theoretical tool for the whole of mode coupling and conversion phenomena 
(except a few cases when rigorous analytical treatment is possible j^) is the coupled-mode, 
or the coupled-wave, formalism 3, 0] that was first proposed more than fifty years ago and 
since then became commonplace in both optical engineering research and textbook litera- 
ture. The mathematical foundation of all the modifications of the coupled-mode formalism 
is the method of variation of independent coefficients — one of the fruitful methods in 
perturbation theory 0. Nevertheless, the formalism grounds upon the differing amount 
of approximations — e.g., that of slowly-varying dielectric function e(r), so that the elec- 
tric field E(r) is assumed to satisfy V(V ■ E) = 0; the so-called parabolic approximation; 
and others — which certainly limits the usefulness of the model. Moreover, the resultant 
"coupled modes" often fail to satisfy orthogonality relations and boundary conditions in the 



2 



actual (perturbed) structure. Finally, the coupled-mode formalism ignores all of the ex- 
plicit mathematical parallels to the classical quantum-mechanics problems that also exploit 
perturbation techniques, in particular the widely known perturbation method for station- 
ary levels 0, — a modification of the Lindstedt-Poincare method of strained parameters 
sometimes called the Rayleigh-Schrodinger method (see Refs. jsll?! for the bibliographical 
details). Meanwhile, establishing and tracking the analogies between classical quantum me- 
chanics and optics tends to be a really stimulating approach — both through adopting the 
formalism of electromagnetic-wave propagation in dielectric materials to coUisionless propa- 
gation of effective-mass electron waves in semiconductor crystals 0], and on the other hand 
through applying the formalism of quantum-mechanical solid-state electronics to electro- 
magnetic propagation in periodic dielectric structures which proved to be a very successful 

,3- 



programme during the last decade 
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In this article, I translate the classical stationary perturbation technique — i.e., the 
method of strained parameters — to the case of mode coupling and conversion in optical 
waveguiding systems. The problem of source-free light propagation is formulated in lSec. lH 
as an ordinary eigenvalue problem in the squared free-space wavenumber k = k"^ = lo^ j(? 
(or, actually, in the angular frequency lS) for the magnetic field H(r), and perturbations 
in dielectric index profile of the waveguide are considered, embracing the cases of optical 
inhomogeneity, magnetization of arbitrary orientation, and nonlinearity. Both the eigen- 
values and the eigenvectors are expanded into the series in perturbation parameter, giving 
rise to an electromagnetic counterpart of the quantum-mechanics stationary perturbation 
method as an alternative to the coupled-mode formalism. Although physical interpretation 
of the former in the case of nondegenerate spectrum of eigenvalues looks rather questionable 
from the point of experiment, the method however provides promptly solvable and naturally 
interpretable treatment of the modal behavior at the anticrossings, which is an issue in the 
mode conversion analysis. To demonstrate the efficiency of the method (to distinguish from 
the coupled-mode theory, the qualifying term "stationary" is retained here with the primary 
meaning that the coefficients of an expansion into which an input wave is decomposed are 
not varying with coordinates), I refer in lSec. IlTl to the two practice-targeted examples from 
integrated optics concerning TE ^ TM mode coupling: a multilayer gyrotropic waveguide 
subject to constant, arbitrarily oriented magnetic field, and an optically-inhomogeneous 
planar dielectric waveguide — an example which falls out of the scope of the standard 
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coulped-mode theory. The unperturbed basis for the both cases is a multilayer waveguide 
composed of linear isotropic dielectric materials — a simple structure known to exhibit under 
certain circumstances the perfect phase matching of guided TE and TM modes thus it 
is natural to implement here the developed theory for the mode conversion at anticrossings. 
For both systems, I pursue the formulated way to derive closed-form first-order analytic 
expressions for the minimum value of mode mismatch and hence for the length of complete 
TEn <-> TMm mode conversion (the coupling length, or the beat length). Finally, I conclude 
in ISec. IVI with some claims concerning further possible applications and modifications of 
the stationary perturbation technique for the optical waveguiding theory. 

II. H-EIGENPROBLEM AND ANTICROSSINGS 
A. H-eigenproblem 

Let me start with the wave equation for the magnetic field 

V X (e^^V X H) = kU, (1) 

where e = e(r) is the dielectric permittivity distribution, k = k"^ = uP' jc? the squared free- 
space wave number, uj the angular frequency, c the vacuum speed of light, and the field H 
satisfies additionally 

V ■ H = 0. (2) 

Eq. (Q) can be treated as an ordinary eigenvalue problem in k (or, actually, in a;) for the 
field H, with Maxwellian operator M. defined by 

A4H = V X (e-^V X H) . (3) 

In fact, a similar equation can be written for, say, the E field: 

e-^V X (V X E) = kE, (4) 

but there are two sound reasons to restrict oneself to the H-eigenproblem (^, but not to 
its E-counterpart (Q: first, an accompanying to Eq. divergence equation (j2I) is clearly 
simpler than V ■ (eE) =0 — a satellite equation for the E-eigenproblem; second, given the 
dielectric index a real scalar function of coordinates — i.e., a real scalar field — or, more 
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generally, a Hermitian dyadic field, the Maxwellian Q becomes Hermitian too |l3|] — the 
fact that is though not crucial in macroscopic electrodynamics (non-Hermitian Maxwellian 
would generete a set of nonorthogonal eigenmodes and complex-valued eigenfrequencies — a 
situation normally faced in optics of lossy qyrotropic media Q|) nor pertinent to the present 
treatment, but aesthetically pleasant and legitimates many of the cross-references between 
Eq. and the stationary Schrodinger eigenproblem with Hermitian Hamiltonian as well. 

For a large class of waveguiding systems, namely, for those exhibiting continuous transla- 
tion symmetry along the direction of light propagation (say, the z direction) we can assign 
harmonic dependence of the fields on that direction; in particular, the magnetic field reads 

H(r;/3) = e'^^?i(x,y), (5) 

where /5 is the propagation constant, so that for a given z-independent dielectric index profile 
e(o)(a;, Eq. (jT)) reduces to an eigenvalue problem for 7i{x,y); the reduced eigenproblem 
operator depends parametrically on (3 then and at a fixed (3 spawns a set of eigenvalues 
thus yielding dispersion structure V{k,[3) = of a perfect waveguide. If an actual 
dielectric permittivity distribution e(r) of the waveguide happens to differ from e(o)(x,?/) for 
which Eq. © holds, that is e(r) = e(o)(x,?/) + £:5e(r) (we do not need to specify at this stage 
whether (5e(r) is a scalar or a tensor field) and 

e-^ (r) = ^(0) (x, y) + £^(1) (r) + £^^(2) (r) + . . . , (6) 

where ^(i)(r) = [— 5e(r)]Y[e(o)(2;, y)]*+^, i = 0,1,..., then in the spirit of the Lindstedt- 
Poincare perturbation method of strained parameters 0, ch. 3] one can rewrite Eq. 
as 

(X(0)+^_^(i)+^2^(2) + , , )H = fi:H, (7) 
A^«H = V X (e(i)V X H), (8) 

in terms of the unperturbed Maxwellian Ai^^^ and the higher-order perturbation operators, 
and expand then the eigenvectors H„ of the perturbed problem into a series in e: 

H„ = HW+5H«+52H(2) + 0(e^), (9) 

and similarly for the eigenvalues: 

^^ = ^m+sK^^^+e\^^^+0{e'). (10) 



The procedure for further solving Eq. ((Tj) using expansions © and (fTUj). i.e., for finding the 
unknown functions H^*'' and is quite well-developed in perturbation theory. Substituting 
Eqs. (fTUI) into Eq. ((Tj) yields for the zeroth-order and linear in e terms: 

^(o)H(o) = 40)H(0), (11) 

^(°)H« + ^(i)Hf = 4°)hW + «:«Hf , (12) 

with the homogeneous boundary conditions V ■ Yl^n = 0. The solution to the unperturbed 
problem pi|) is assumed to be known; to say more — and this is an important point — 
I further assume all the H^^^ vectors to be of the form (0) which is likely to embrace all 
the cases of practical interest. The subsequent steps of the method depend essentially on 
whether the spectrum of AA^'^'^ is degenerate or not. 



B. Nondegenerate spectrum 

If there are no degenerate eigenvalues among those of A^'-'^-' spectrum, then we customarily 
obtain up to the first-order correcting terms y]: 



H„(r; /3) ^ HW(r; /3) + . fi^HW(r; /5) (13) 



and 

«„(/?) ^4°)(/3)+£A<«(/3), (14) 
where the first-order coupling matrix 

Mi^li(3)^M-'[ U<'\v;P)M^'\r)U^^\r;(3)dv. (15) 
Jw 

The integration volume W embracing the waveguide is formally infinite, dv = dxdydz, 
the normalizing factor Af equals either Z — the distance of light propagation within a 
waveguide in the z direction (i.e., the length of the waveguide) — in the case of cylindrical 
(two-dimensional) geometry, or S* = YZ — the yz square of a waveguiding cell — for planar 
(one-dimensional) systems. This normalization is due to specific orthogonality condition 
assigned to the Hn ^ modes: 

w Jw 

5nmZ J \Hm{x, y)f dx dy = SnmZ (cylindrical waveguide), 

(16) 

SnmS J \7im{x)\ dx = 6nmS (planar waveguidc). 
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Normalized wavenumber kK 

FIG. 1: Modal dispersion curves ( ) of the unperturbed hypothetical structure and /? = const 

lines ( ) for the case of (a) nondegenerate spectrum of M.'^^^ and (b) the two modes, H^''^ 

and , exhibiting degeneration. Here A is a characteristic dimension of the structure; the grey 
circles mark the modes related to the expansions and (|17|) . 

I explicitly designated the (parametric) dependence of all the quantities in Eqs. (O, 
(HD), and (fT3j) on [3 entailed by the previous assumption that the unperturbed modes have 
the form as per Eq. (jSj). Now, the restriction (3 = fixed for the set of nondegenerate modes 
Hi°^e~*'^"* oc e*^^^"'^"*) = Q^^r^i^nz-ct) jj^^Q which the nth perturbed eigenmode H„ is expanded 
via Eq. — modes characterized thus by essentially different eigenwavenumbers kn and 
differing normalized propagation constants Pn = P/kn — seems to correspond to an utterly 
odd situation from the point of experiment (see |Fig.~H ). To aid this shortcoming somehow, 
one might consider Eq. (0) as a generalized eigenproblem in, e.g., (3 or j3, parametrically 
dependent on k, which would lead to expansions like (0 but into the {Hi°^(r; k)} set. Of 
course, this severely complicates the formalism — and probably here is an issue why the 
coupled-mode theory has been exclusively dominating for decades in the field. Fortunately, 
the case of quasi-degenerate eigenvalues allows promptly solvable and easily interpretable 
treatment in terms of the ordinary H-eigenproblem in k. 
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C. Degenerate spectrum and anticrossings 



Let us focus on the case of two-degenerate or nearly degenerate eigenmodes H^°'' and 
Hg^''. Borrowing the known result from quantum mechanics j^, we see that in this case 
both the H^°'* and 112''^ modes dominate the expansion (jU)) of the related fundamental modes 

H± = Ctnf'^ + Ca^Hf + 0{e) (17) 

with the constant (in particular, ^-independent) coefficients 



± 



M 
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11 - M22y + 4:Ml2M21 

Mil - M22 



(18) 



(19) 



^/{Mii - M22Y + 4A<i2A^2i _ 
Here and below in the article I pursue only the first-order approximations and hence omit the 
overly superscript ^^'> over the coupling matrix elements A1„m, now /3-independeiitj defined 
by Eq. ()15|). The first-order correction to the corresponding eigenvalues leads to 



4°^ + 4°^ + e{Mii + M 



22 



± 



K 



(0) 



K 



M 



22 



+ e^Mi2M2i. (20) 



In quantum mechanics, this situation corresponds to quasi-degenerate energy levels of a 
quantum system and is associated with enormous variety of effects and treatments — take 
the Landau-Zener model for example |15|, nourishing nearly 100 publications a year. In 
the realm of optics, Eq. (fTTj) with coefficients given by Eqs. (fTHj) . (fT^ is fully consistent 
with the general statement that complete mode conversion can be achieved in any system 
exhibiting anticrossings. Indeed, instead of the h[°'' mode and the H2'^'* mode crossing one 



another on the dispersion structure diagram (Fig. 1 ), at the anticrossings the H^^"* mode 



"transforms" continuously into the H2°^ one, and vice-versa, over a certain frequency range. 



so it is intuitive to expect that somewhere in the vicinity of the virtual crossing point the 
H^'^'' and Hg^"* modes would be equally unsuitable to play the role of the fundamental modes 
H+ and H_, and if the original (input) wave is Hf^ or Uf, it wouldn't "feel" which of the 



fundamental modes, H+ or H_, is more suited for it, and thus would oscillate between 
and ^2^ . Mathematically, these heuristic speculations are expressed via the requirement 

1 



(0) 



ic+i = |cr 



(21) 
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(under appropriate normalization) necessary for the complete mode conversion to occur, 
which is equivalent, as immediately follows from the comparison of Eqs. ()18|) and (fTI?|) . to 
canceling out the A4u and elements under the root in Eqs. (fTH|) . (fT^ . and due to 

Mn = M22- (22) 

This condition does not necessarily correspond to the case of stringently degenerate eigen- 
values and h!"2 \ but with a good deal of reason we may assume h!^^ — k^2^ = 0(£:^) to 
estimate the minimum mode mismatch Akmin from Eq. ()2()|1 . With this assumption, it is 
under condition that the quantity 

where ko = ujq/c, uq is the frequency of a virtual mode crossing (phase matching), reaches 



its minimum value of 2ey/jAi2M^, thus the minimum eigenmode mismatch Akmin along 
the f3 = const curve corresponds to the complete mode conversion and reads 

= €koWMi2M2i. (24) 

In a nonabsorbing waveguide (in Sec. IIIII we will deal with the two) Ainm is a Hermitian 
matrix, hence = -^21 ^'i- is simplified to 

Akmin = ekQ^\Mi2\. (25) 

Complete mode conversion occurs when the and modes propagating in the z 
direction accumulate a phase difference of vr: AkminP^c = vr; substituting here Eq. ((23), we 
obtain for the corresponding coupling length: 

ef3\Mi2\ 

where, as before [see Eq. fll5|) ]. the coupling matrix element 

M,2 = M-^l [V X (e(i)V X Hf))] ■ YifUv. (27) 
Jw 

In contrast to the coupled-mode theory, the M.12 element might and, in general, does depend 
on the length of a waveguide along the z direction because of the possible z dependence of 
the (^(1) function. In what follows I calculate M.12 and hence Zc explicitly for the two specific 
systems relevant to integrated optics: for a multilayer gyrotropic stack and for an optically- 
inhomogeneous waveguide. 



III. MODE CONVERSION AT PLANAR GEOMETRY 



A. Gyrotropic waveguide 

Let the unperturbed waveguide be described by the scalar dielectric permittivity dis- 
tribufou As argued by symmetry reasoning, relafvity eo,rsideratio,rs. and energy 

conservation [13|, the most general form to allow for the electromagnetic disturbance 
induced by the eh field applied to (initially) isotropic medium is 

e^ = e(o)T + i£Ch'' +£'r/(h®h-h2), (28) 

where I is a unit three-dimensional dyadic; ( and t] are real functions of h^; the superscript 
^ over the h vector denotes the antisymmetric dyadic dual to h, so that h^g = h x g gives 
conventional vector product of h and g; the circled cross ® denotes the outer product — 
the dyad: h«g = k.,„ = 1,2,3. If . . small, we may assmrre C and , to be constant 
(this assumption leads to the introduction of gyration vector [13|, |1J]) and write down for 
the inverse to the e* dyadic jl4| 

^-i = e-;T-^£u><+0(£2), (29) 
in terms of u = e^^^h, so that the coupling matrix element can be written as 

Mil' = -kS-' I [V X (u X Df ))] ■ Ylf^dv, (30) 
Jw 

where I used the frequency- domain Maxwell's equation for V x H in the source-free region, 
V X H = — i/cD. Equation (jHUj) can be significantly simplified by applying the vector 
identity [3| 

V X (u X D) = u(V ■ D) - D(V ■ u) + (D ■ V)u - (u ■ V)D (31) 

valid for any three- vectors u and D. The first term on the right vanishes immediately due 
to V ■ D = 0; then, in a specific but physically reasonable and normally considered case of 
u being a constant (or piecewise-constant) vector in a given volume V and equaling zero 
outside of it, the second and third terms of the sum in Eq. being multiplied by H and 
integrated as per Eq. (jHIH) lead to the surface integrals instead of the volume ones. For the 
second term one obtains: V ■ u = — (fi ■ u) 5{r — R), where R = R(r) defines the boundary 
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Ti of V, h is the outer normal to S; hence 

/ [D( V ■u)]-Hdv = - I [D(n ■ u)] ■ H da. (32) 
For the third term we similarly have: (D • V)u = — (D ■ n)u5(r — R), hence for an integral 

f [(D • V)u] ■Yldv = - [ [(D • n)u] ■ H da. (33) 

Jw Jt, 

Thus for the case of constantly magnetized waveguide (up to this point, no assumptions have 
been made regarding its actual geometry) the integral in Eq. (jHUjl is reduced via Eqs. (jSH), 
dSl, and dSl to 

[ [Vx (uxD]-Hrft; = [ [D(n-u)]-H(ia- / [(D-n)u] -Hda- /" [(u- V)D] ■ H rfi;. (34) 
Jw Jt, Jt. Jv 

Now let us switch to the planar structures. In a planar unperturbed waveguide the 

eigenmodes are classified in terms of TE, or 5-polarized, and TM, or p-polarized waves (for 

an elegant derivation of this common fact via symmetry considerations see Ref. Q]). In a 

coordinate system in which the x axes is normal to the bimedium interfaces and the light 

energy is guided in the z direction, we can write down, to evaluate the M.^^^ element as per 

Eq. dnni): 

Df = D, = e^^^0,(a;)y (35) 



for the s-polarized mode, and 



Hf = Hp = e^'^^^,(x)y (36) 



for the p-polarized one. Here (j)s{x) and ippix) are the lateral distributions of and Hp 
respectively, obtained from the unperturbed Maxwell's equations for the given virtually- 
phase-matched modes. We see immediately from Eq. ()35|) that the second surface integral 
in Eq. ()34|) vanishes at planar geometry, since ■ n = ■ x = ■ (— x) = holds 
identically in this case. Finally we arrive at 

= -kS-^f [D,{h-u)]-li;da + kS-^ [ [{u-V)JD,]-li;dv 

Jy, Jv 



-k 



J Y \ dx dx J J Y \ d^ 



(37) 



11 



where is the value of vector u in the ath layer bounded between the x = Xa-i and x = Xa 
planes, and the integrals 

j.f<tsd^\ ^ r M^)dM-) (38) 

J-, fm , r *-(f['(-> (39) 

a = 1, . . . , n (ri is the number of layers). 

Eq. (jHTj) accounts for no mode conversion at transversal magnetization, in agreement with 
the long ago established result |13]. A few more implicit tips could also be deduced thereof: 

1. In a symmetrically sandwiched stack, the polar magnetization (u = m^x) virtually 
doesn't couple the modes exhibiting similar symmetry in the and ipp{x) functions 
(that is, when both (j)s{x) and ipp{x) are either even or odd); on the contrary, the modes 
of opposite symmetry in 0s(x) and ipp{x) are virtually not sensitive to the longitudinal 
magnetization (u = u^z). Asymmetric sandwiching impairs this behavior. 

2. The integral (jH^ contains scaling parameter f3~^ that defines characteristic thickness 
of the waveguide corresponding to the comparable values of the both integrals, Eqs. ()38p and 

in Eq. ^7\i. and hence to the comparable shares of polar and longitudinal polarizations 
in the mode conversion efficiency. For optical frequencies f3~^ ~ 100 nm; if the thickness of 
the waveguide considerably exceeds than 'ipp{x)would appear to be a too slowly varying 
function and therefore 

J'(^-^J-(^- (40) 



dx J \ /3 

For the ultrathin layers the inverse inequality holds, but in this regime at most one mode is 
guided in the structure, which is apparently out of the scope here. 



B. Optically-inhomogeneous dielectric waveguide 

If we assume = .^(r) in Eq. (j27j) to be a scalar function within the waveguide volume 
V, which conveys the cases of optical inhomogeneity and "isotropic" nonlinearity, then it is 
advantageous to simplify Eq. (j^Tj) using the identity 

V X (eV X H) = V(Ve ■ H) - (H ■ V) - (V^ ■ V)H - ^V^H, (41) 

where V • H = is taken into account. In planar geometry, we are concerned as before with 
the TE and TM modes, but now expressed exclusively through the magnetic field vectors. 
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For the TM (p-polarized) mode we can readily apply Eq. (j^ : for the TE (s-polarized) mode 
we have for the electric field, in parallel with Eq. ((23), 



E, = e^^>,(x)y, (42) 

hence 

H, = -ik-'V X E, = k-'e'^' (^M(x)x - ^^^^^z) , (43) 

where we once again encounter /3 (in fact, in the role of scaling parameter. If we put, 
say, H = Hs in Eq. (PT|) and multiply the result by , then only the first and second terms 
of the sum (j4H) will produce nonzero results, as it immediately follows from Eqs. (j43p and 
flHUI) . Equation ^T7\i thus gives 

M:f = S-^! [V(V^H,)]■H;rf^;-5-l/ [(H,- V)V^]-H;rft;. (44) 
Jw Jw 

Now care should be taken when dealing with V.^ function which is discontinuous in the 

vicinity of material boundaries [otherwise — if is assumed to be all-continuous — Eq. 



gives M.'^^ = 0]. The safest way is to introduce a new vector field <; = <j(r) instead of V^, 
with a requirement V x <j = everywhere except at the boundaries: 

V X = -(n X <j)5(r-R), (45) 

and with this in mind to perform the brute-force evaluation of integrals in Eq. ()44|1 : 

/ [ViV^-U,)]-U;dv = pk-'[ ^^^{x)ij,{x) dv - zk^' [ ^^^ij^{x)dv, (46) 

Jw Jw Jw '^y 

/ [(H,- V)Ve]-H;dt; = /3fc-i/ ^^^(x)ij,{x)dv-zk-'[ ^^^^p^(x)dv, (47) 
Jw Jw '^•^ Jw ^-^ 

hence 

MTJ' = PikSy [ {-Vx^),^,{x)^j,{x)dv + z{kS)-'[ {-V x<;),^^^^p{x)dv 
Jw Jw 

= (3{kSy^ 'iy{r)(ps{x)ipp{x) da 



s 

,-11 I c'^i^fi^^PJ , c» 



V 



ms)-' I { e;^^^^ + ) dv. (48) 



We see that for the nonzero mode conversion, the ^(r) function should exhibit explicit 



dependence on the y coordinate. For this reason would depend on the transversal 



sp 



dimension Y of the waveguide and on the values of .^(r) at the y = and y = Y boundaries. 
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Another important consequence relates the effect of the z modulation of ^ (r) on the coupling 
length. If ,^(r) does not depend on z, then Ailp^ would in turn be independent on the waveg- 
uide length Z; but on the contrary, whenever ^(r) is a stochastically oscillating or periodic 
function of z, the coupling matrix Ai^^ and hence the coupling length Zc would become 
Z-dependent. Say, for the ^{z) oc co?,{z/zq) dependence we have J^'^^{Z) oc Z^^smi^Zlzo) 
which is of the order of Z~^ and brings about an unexpected scaling rule z^ ^ Z for the 
coupling length in a planar waveguide with periodic modulation of 5e(r) in the z direction. 

Finally, I would like to remind here that conventional expression for the off-diagonal 
coupling matrix element given by the standard coupled- mode theory 



gives an identical zero for the TE ^ TM mode conversion, being the dielectric perturbation 
e 5e(r) a scalar function. 

IV. CONCLUSION 

Mode conversion displays breaking the (initial) symmetry of the Hamiltonian of a system 
by some perturbation that distorts the mode spectrum, i.e., shifts the eigenvalues and alters 
polarization of the eigenmodes, hence a natural tool to treat the mode coupling and conver- 
sion phenomena is the perturbation technique. For the mode conversion at anticrossings, I 
presented in this article an electromagnetic counterpart of quantum-mechanical perturba- 
tion theory for quasi-degenerate levels based on the Lindstedt-Poincare method of strained 
parameters as a specific alternative to the entrenched coupled-mode formalism grounded 
upon the method of variation of independent coefficients. 

The general expressions derived for the minimum mode mismatch at anticrossings, 
Eq. (j2ni), and for the coupling length, Eq. ^I^. are compact, transparent and premise on the 
calculation of just one element of the coupling matrix, Eq. ()27|1 . That matrix element was 
calculated explicitly for the two cases of interest in integrated optics: for a multilayer gy- 
rotropic waveguide under piecewise-constant, arbitrarily oriented magnetization [Eq. (|H7jl ]. 
and for an optically- inhomogeneous planar dielectric waveguide [Eq. PHj) ]. In a similar way 
a large variety of optical systems can be analyzed, including Bragg's fibers, photonic crystals 
with broken periodicity, etc. 




(49) 
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Finally, I should note that by means of an appropriate coordinate mapping, the problem 
of perturbation due to shifted material boundaries — so to say, geometrical perturbation 
known to spur difficulties when treated via conventional perturbation techniques — 
can be reduced to the problem of perturbation in permittivity and permeability profiles of 
a waveguide exhibiting perfect geometry in those (curvilinear) coordinates and governed 
by "Cartesian-looking" Maxwell's equations. Such a trick that recently aided the FDTD 
modeling of high index-contrast photonic crystals [19] would also significantly augment both 
the standard coupled-mode theory and the stationary perturbation technique described here. 

Acknowledgments 

I cordially thank Prof. I. D. Feranchuk who gave an initial push for this work by making 
me look at the whole variety of quantum-mechanical phenomena associated with quasi- 
degenerate levels. 



[1] D. Marcuse, Theory of Dielectric Optical Waveguides (Academic Press, New York, 1974); 

A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapmen&Hall, London, 1983); 

A. Yariv and P. Yeh, Optical Waves in Crystals: Propagation and Control of Laser Radiation 

(Wiley, New York, 1984), and references therein. 
[2] T. D. Engeness, M. Ibanescu, S. G. Johnson, O. Weisberg, M. Skorobogatiy, S. Jacobs, and 

Y. Fink, Optics Express 11, 1175 (2003). 
[3] S. G. Johnson and J. D. Joannopoulos, Photonic Crystals: The Road from Theory to Practice 

(Kluwer, Boston, 2002). 

[4] S. Zhang and T. Tamir, Opt. Lett. 20, 803 (1995); S. Zhang and T. Tamir, J. Opt. Soc. Am. 

A 13, 2403 (1996). 
[5] W.-P. Huang, J. Opt. Soc. Am. A 11, 963 (1994). 
[6] A. H. Nayfeh, Perturbation Methods (Wiley, New York, 1973). 

[7] C. H. Wilcox, Perturbation Theory and its Applications in Quantum Mechanics (Wiley, New 
York, 1966). 

[8] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory (Pergamon 



15 



Press, New York, 1965), 2nd ed. 
[9] G. N. Henderson, T. K. Gaylord, and E. N. Glytsis, Phys. Rev. B 45, 8404 (1992); T. K. 

Gaylord, G. N. Henderson, and E. N. Glytsis, J. Opt. Soc. Am. A 10, 333 (1993). 
[10] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987); S. John, Phys. Rev. Lett. 58, 2486 (1987). 
[11] J. D. Joannopoulos, R. D. Meade, and J. N. Wmn, Photonic Crystals: Molding the Flow of 

Light (Princeton University Press, Princeton, NJ, 1995). 
[12] D. M. Shyroki and A. V. Lavrinenko, J. Opt. A: Pure Appl. Opt. 5, 192 (2003). 
[13] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media (Pergamon Press, 

Oxford, 1960). 

[14] F. I. Fedorov, Theory of Gyrotropy (Nauka i Technika, Minsk, 1976). 
[15] V. A. Benderskii, E. V. Vetoshkin, and E. I. Kats, farXiv : cond-mat/0303275| (2003). 
[16] G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers (McGraw- 
Hill, New York, 1961). 
[17] J. R. Gillies and P. Hlawiczka, J. Phys. D: Appl. Phys. 9, 1315 (1976). 

[18] N. R. Hih, Phys. Rev. B 24, 7112 (1981); S. G. Johnson, M. Ibanescu, M. Skorobogatiy, 
O. Weisberg, J. D. Joannopoulos, and Y. Fink, Phys. Rev. E 65, 066611 (2002). 

[19] A. J. Ward and J. B. Pendry, J. Modern Opt. 87, 8268 (1996); Phys. Rev. B 58, 7252 (1998); 
Comput. Phys. Commun. 128, 590 (2000); D. M. Shyroki, .arXiv : physics/0307029 (2003) . 



16 



